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Abstract 

Background: Malaria remains a serious problem in French Guiana, which is at potential risk for drought linked with 
the El Nino Event and where there could be a risk of malaria epidemic after the onset of an El Nino event. 

Methods: A time series analysis using ARIMA was developed to investigate temporal correlations between the 
monthly Plosmodium falciparum case numbers and El Nino Southern Oscillation (ENSO) as measured by the 
Southern Oscillation Index (SOI) at the Cayenne General Hospital between 1996 and 2009. 

Results: The data showed a positive influence of El Nino at a lag of three months on P. falciparum cases (p < 0.001). 
The incorporation of SOI data in the ARIMA model reduced the AIC by 4%. 

Conclusions: Although there is a statistical link, the predictive value of ENSO to modulate prevention intervention 
seems marginal in French Guiana. However, additional work should refine the regional dependence of malaria on 
the ENSO state. 



Background 

In French Guiana, malaria remains a serious problem 
with a number of malaria cases observed every year ran- 
ging from 3,500 to 4,500. Both Plasmodium falciparum 
and Plasmodium vivax are prevalent. Since 2000, 60% of 
the cases are due to P. vivax [1]. 

The ENSO phenomenon refers to the cyclic warming 
and cooling of the equatorial Pacific Ocean coupled 
with changes in the atmospheric pressure across the 
Pacific. This is the most important climatic cycle contri- 
buting to worldwide interannual variability in climate 
and the likelihood of climatic anomalies. The two 
extremes of ENSO are El Nino (a warm event) and La 
Nina (a cold event), which can create rainfall and tem- 
perature fluctuations [2]. 

The El Nino-Southern Oscillation (ENSO) is the best- 
known example of quasi-periodic natural climate 
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variability on the interannual time scale. It comprises 
changes in sea temperature in the Pacific Ocean (El 
Nino) and changes in atmospheric pressure across the 
Pacific Basin (the Southern Oscillation), together with 
resultant effects on world weather. El Nino events occur 
at intervals of 2-7 years. In certain countries around the 
Pacific and beyond, El Nino is associated with extreme 
weather conditions that can cause floods and drought. 
The other extreme (La Nina) (cooling) is the reverse 
pattern of climate anomalies. 

There is strong evidence that El Nino Southern Oscil- 
lation (ENSO) is associated with increased risk of 
malaria in regions of the world where climate is linked 
to the ENSO cycle [3]. These include, among others, 
countries in South America [4-8]. At a regional level, 
French Guiana is located in a potential risk area for 
drought linked with the El Nino Event and where there 
could be a risk of malaria epidemic after the onset of an 
El Nino event [3]. 

The purpose of the investigations reported here was to 
test and quantify the influence of ENSO on malaria 
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transmission in French Guiana using monthly data of 
malaria cases at the Cayenne General Hospital Because 
of the difficulty to discern P. vivax malaria relapses 
from re-infections, this study focused only on P. falci- 
parum malaria transmission. 

Methods 

ENSO data 

The Southern Oscillation Index (SOI) is designed to 
measure the strength and phase of the Southern Oscilla- 
tion. SOI is a measure based on the differences in the 
atmospheric pressure between Tahiti in the eastern equa- 
torial pacific and Darwin in Australia (West Pacific), 
expressed as a standard deviation form the norm. SOI is 
used to quantify the strength of an ENSO event. It is gen- 
erally negative during an El Nino event and positive dur- 
ing a La Nina event. SOI was taken from the Australian 
bureau of meteorology [9] . 

Malaria data 

In spite of some encouraging results to statistically dif- 
ferentiate P. vivax malaria relapses from re-infections at 
individual scale [10], due to the difficulties to discern 
P. vivax malaria relapses from re-infection in endemic 
areas as French Guiana [11], this study focused only on 
P. falciparum transmission. 

All clinical malaria episodes diagnosed for patients 
consulting the emergency service of Cayenne hospital 
were retrospectively collected. Malaria was defined as 
temperature > 38°C at the time of consultation or fever 
within the past 48 hours associated with Plasmodium 
asexual forms on a thin blood smear (screening sensitiv- 
ity: 50 Plasmodium per \A) 

Data were extracted from the Parasitological Labora- 
tory files using designated software in accordance with 
legal procedures concerning confidentiality. Malaria data 
from 1996 to 2009 were recorded as numbers of cases 
of P. falciparum malaria per month. 

Data analysis 

ARIMA models were used to evaluate the relationship 
between SOI and monthly malaria cases. Before, the 
ARIMA methodology application, a log-transformation 
of the series was made to stabilize variance. 

Modelling with ARIMA involves the estimation of a 
series of parameters to account for the inherent dynamics 
in the time series, including the trends and autoregressive 
and moving average processes. The general model intro- 
duced by Box and Jenkins [12] includes autoregressive 
and moving average parameters, and explicitly includes 
differencing in the formulation of the model. An ARIMA 
(p, d, q) model comprises three types of parameters: the 
autoregressive parameters (p), number of differencing 
passes (d), and moving average parameters (q). The 



multiplicative seasonal ARIMA (p, d, q) (P, D, Q)s model 
is an extension of the ARIMA method to time series in 
which a pattern repeats seasonally over time. Analogous 
to the simple ARIMA parameters, the seasonal para- 
meters are: seasonal autoregressive (P), seasonal differen- 
cing (D), and seasonal moving average parameters (Q). 
The length of the seasonal period is represented by s. 

Various permutations of the order (p, d, q) (P, D, Q)s 
were computed, and the optimal combination of para- 
meters using the « Akaike's information criterion (AIC) 
» was chosen. The smallest values of Akaike's informa- 
tion criterion (AIC) were set as the standard to identify 
the best-fit model [12]. The correlogram and partial cor- 
relogram graphs were also used to help in fixing the 
value of orders to include in the model. The likelihood 
ratio test was used to determine if inclusion of other 
covariates helped to improve the fit of the model. 

In the model, SOI, at lags of 0 to 12 months, respec- 
tively, was then fitted into the ARIMA model to screen 
for potential predictors of the P. falciparum case num- 
bers. Lag series significantly associated with malaria, 
(p-value <0.05), were singled out to fit the best multivari- 
ate ARIMA model. The Ljung-Box Q test was applied to 
ascertain whether the residual series were white noise. 

Results 

The monthly observed cases of P. falciparum malaria, 
and SOI, were shown in Figure 1. Of all the models 
tested, the ARIMA (1, 1, 1) model fitted the data best 
according to AIC. The best fitting model showed that no 
significant seasonal pattern explained the monthly cases 
of P. falciparum malaria alone. The results of ARIMA 
analysis reveal that the logarithm of the malaria cases 
was significantly influenced by SOI (Table 1). 

In this model, the local moving average parameter is 
-0.8961 (p =<0.001), and the auto-regression term is 
0.3813 (p =<0.001). With SOI data in the model, loga- 
rithm of P. falciparum cases number was negatively asso- 
ciated with SOI at a lag of three months (P = -0.0121, p = 
0.006). The incorporation of SOI data in the ARIMA 
model reduced the AIC by 4%. 

Discussion 

There were two main limitations in the present study. 
The first one is that this study has used aggregated data 
coming from the whole heterogenic territory of French 
Guiana and has to be completed by additional studies at 
a local scale to determine whether a greater proportion 
of the variance is predicted, and thus giving a more 
practical value to the observed correlations in terms of 
prevention. 

The second one is that, although 60% of malaria cases 
are caused by P. vivax in French Guiana, P. vivax cases 
were not included in this analysis due to impossibility to 



Hanf et al. Malaria Journal 201 1, 10:100 
http://www.malariajournal.eom/content/10/1/100 



Page 3 of 4 



Evolution of P. falciparum cases 




1996 1998 2000 2002 2004 2006 2008 2010 
Year 



Evolution of the Southern Oscillation Index (SOI) 




— i 1 1 1 1 1 1 1 — 

1996 1998 2000 2002 2004 2006 2008 2010 

Year 



Figure 1 Plasmodium falciparum cases at the Cayenne General 
Hospital plotted with the Southern Oscillation Index (SOI), 
1996-2009. 
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differentiate relapses from re-infections due to missing 
detailed patient follow-up data. 

Nevertheless, this study showed a consistent temporal 
correlation of ENSO with P. falciparum cases in Cayenne 
General Hospital with a three-month time lag from 1996 
through 2009. Although ENSO influence on malaria 
transmission was suspected in French Guiana, this study 
is the first to identify a statistically significant correlation 
between El Nino and malaria cases in this area. A 
decreasing in SOI three months before was found to act 
as a positive factor on P. falciparum transmission. 

The identification of a three-month lag allowed to 
quantify the delay between ENSO perturbation and con- 
sequences on malaria epidemics in French Guiana. 



Table 1 Estimated coefficients by ARIMA Model between 
SOI and logarithmic transformation of malaria cases in 
Cayenne General Hospital, 1996-2009* 





Coef. 


Std. Err. 


P.value 


AR 


0.3813 


0.0941 


< 0.001 


MA 


-0.8961 


0.0479 


< 0.001 


SOI lag 3 


-0.0121 


0.0049 


< 0.001 



AR: Autoregressive parameter 
MA: moving average parameter 
AIC of the model: 261.12 



Possible explanations for the identified association 
between ENSO and malaria in French Guiana include 
the effect of climate on the population dynamics of vec- 
tors (mainly Anopheles darlingi in French Guiana) 
through changes in population densities or survival 
rates, through availability of adequate breeding sites, but 
also through parasite associated parameters and man- 
vectors contacts. 

The data indicated that ENSO explain 4% of the varia- 
tion in malaria case numbers in the studied region. At 
the study scale, the remaining 96% of the variation pos- 
sibly involved non-climatic causes, such as population 
immunity and socio-environmental factors that influ- 
ence the breeding and ecology of mosquito vectors. This 
result suggested that changes in French Guiana are 
more likely to be a reflection of the efficiency of local 
control measures, population movements, therapeutic 
change, and others social factors [13]. In French Guiana, 
the role of recent immigration for gold-mining activity 
is well known as a main drivers of malaria epidemics 
[14]. 

To conclude, although there is a statistical link, the 
predictive value of ENSO to modulate prevention inter- 
vention seems marginal. However, the weak predictable 
consequence of ENSO for malaria transmission is an 
important piece of information for the development of 
future early warning systems able to detect efficiently 
malaria epidemics in the area. 
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